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ABSTRACT 


Values  of  system  availability  are  computed  for  a  system 
whose  failure  density  is  exponential  (A)  and  whose  repair 
density  is  special  Erlang  (y,a).   The  system  is  modeled  as 
an  alternating  renewal  process.   The  Laplace  transform  of 
the  availability  function  is  developed  and  inverted  (by 
using  a  numerical  procedure)  to  obtain  the  values  of  the 
system  availability. 

Corresponding  values  of  the  long-run  average  system 
availability  and  the  average  availability  over  the  first 
mission  time  are  also  computed.   Comparisons  are  made 
which  establish  that  using  the  long-run  average  value  to 
approximate  the  true  availability  over  the  first  mission 
time  is  a  conservative  practice. 
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I.   INTRODUCTION 

A.  THE  SYSTEM 

Most  physical  systems  people  encounter  are  of  the  type 
which  operate  for  a  time,  fail,  are  repaired  and  put  back 
into  operation.   The  system  then  may  be  said  to  have  two 
states.   Truely,  each  of  these  states  may  be  the  union  of 
numerous  sub-states,  but  the  two  major  states  prevail. 

Some  usual  terms  applied  to  the  first  state  are 
'operational',  'available'  or  'up'.   This  is  the  state  or 
condition  of  the  system  when  it  is  capable  of  accomplishing 
the  task  for  which  it  was  designed  and  built.   The  second 
state,  the  compliment  of  the  first,  is  the  less  desirable 
state  and  is  usually  termed  'failed',  'nonoperational ' , 
'unavailable'  or  'down'. 

Hereinafter  the  system  states  will  be  referred  to  as 
' up '  or  ' down ' . 

B.  DEFINITION  OF  AVAILABILITY 

A  good  deal  of  effort  is  being  expended  by  the  U.S. 
Government  and  by  private  industry  to  obtain  measures  of 
certain  quality  characteristics  applicable  to  such  systems 
as  is  described  in  A.   Hosford  [2]  has  discussed  several 
of  these  characteristics. 

One  such  quality  characteristic  that  is  of  primary 
concern  is  system  availability.   The  'system  availability' 
or  simply  'availability'  is  the  probability  that  the 


system  is  up  at  any  time  t  2.  0 ,  given  it  was  up  at  time 
zero . 

The  availability  function,  n(t),  is  a  real  valued 
function  which  gives  this  probability  for  any  time  t  _>0. 

C.   METHOD  CURRENTLY  PRACTICED  TO  DETERMINE  AVAILABILITY 

In  reliability  literature,  availability  is  often 
defined  differently  as: 


=     MTTF 


MTTF  +  MRT 


where  MTTF  is  Mean-Time-To-Failure  and  MRT  is  Mean-Repair- 
Time.   This  definition  provides  the  only  method  currently 
in  wide  spread  use  for  determining  availability. 

Note:   This  quantity  has  the  variable  name  ZBAR  in  the 
computer  program  and  output  (Appendices  A  and  C). 

Cox  [1]  has  shown  that  regardless  of  the  distributions 
of  uptimes  and  downtimes 


lim  n(t)  =  A 

t  -*■  °° 


That  is,  A  is  actually  the  steady-state  system  availability 

Similarly,  it  may  be  shown  that  the  long-run  average 
availability 

1  t 

lim  -r  f      n  (x)dx  =  A. 

t  +  °°   0 


Jacobson  [3]  describes  a  statistical  procedure  for 
obtaining  a  confidence  interval  estimate  for  A. 

D.   OBJECTIVES  OF  THIS  RESEARCH 

The  objectives  of  this  research  were  the  following: 

1.  To  compute  exact  values  of  the  availability  function 
for  a  system  with  given  distributions  of  uptimes  and 
downtimes . 

2.  To  compare  true  availability  of  a  system  with  that 
obtained  by  methods  currently  in  practice. 


II.   SUMMARY  AND  CONCLUSIONS 

A  physical  system  having  an  'up'  state  and  a  'down' 
state  and  with  the  ability  to  transit  between  the  two 
states  was  modeled  as  an  alternating  renewal  process  for 
the  purpose  of  finding  the  system's  availability  function, 
n(t).   An  explicit  general  expression  for  n(t)  was  not 
found  because  of  the  complex  structure  of  the  model 
equations  under  the  assumed  distributions  of  uptimes  and 
downtimes.   However,  n  (s),  the  Laplace  transform  of  n(t), 
was  developed  from  the  assumed  probability  density  functions 
and  the  renewal  density  functions  of  the  model. 

A  computer  program  was  written  to  use  a  numerical 
procedure  to  invert  the  Laplace  transform  for  specific 
values  of  the  parameters  and  of  time,  the  independent 
variable . 

Upon  consideration  of  the  results  and  some  sample 
problems  (see  Appendices  A  and  B)  the  following  conclusions 
were  made. 

1.  Under  the  fairly  realistic  assumptions  of  the 
model,  exact  values  of  n(t)  can  be  and  were  computed. 

2.  Use  of  long-run  average  availability  to  approximate 
the  true  availability  over  the  first  mission  time,  as  is 
currently  being  done  in  practice,  is  a  conservative 
procedure . 


3.   The  amount  by  which  the  long-run  average  availability 
underestimates  the  true  availability  during  the  first 
mission  time  is  dependent  upon  the  long-run  value,  i.e., 
the  lower  the  long-run  value,  the  more  it  underestimates 
the  true  availability. 


III.   THE  MODEL 

A.   DESCRIPTION 

The  system  was  modeled  as  an  alternating  renewal  process 
with  two  states.   That  is,  the  new  system  was  assumed  to 
start  operating  In  the  up  state.   The  time  spent  in  this 
state  before  the  first  failure,  whether  in  use  or  just 
available  for  use,  was  a  random  variable  LL  .   Upon  failure 
the  system  entered  the  down  state.   The  time  spent  down 
before  the  first  repair,  whether  awaiting  repairs  or 
actually  being  repaired,  was  a  random  variable  D-,  .   The 
length  of  the  first  cycle  was  the  random  time  U,  +  D,  . 

It  was  assumed  the  system  kept  alternating  between  up 
and  down  and  each  repair  returned  it  to  'like  new'  condi- 
tion so  that  the  uptimes  U. ,  i=l,2,...,  were  taken   to  be 
independent,  identically  distributed  random  variables. 
They  were  assumed  to  be  exponentially  distributed  with 
parameter  X   >0  and  probability  density  function 


,  t  <  0 


Note:   All  time  units  in  this  thesis  are  mission  times. 

The  downtimes  D . ,i=l ,2 , . . . ,  were  also  assumed  to  be 
independent,  identically  distributed  random  variables  and 
independent  of  the  uptimes.   They  were  assumed  to  have  a 
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special  Erlang  distribution  with  parameters  y  >  0  and 
01=1,2,...,  and  with  probability  density  function. 


rD(t)  = 


0 


(a-!)! 


r  t 


a-1   -yt 


,   t  <  0 


0  <  t 


The  Special  Erlang  Density  Function 


Weiss  [5],  through  a  simple  argument,  has  shown  the 
renewal  density  functions  for  such  a  renewal  process  to  be 


WQ(t)  =  fy(t)  +  /   W1(t)  fyCt-x)  dx 


W  (t)  =  /   W„(t)  rn(t-x)  dx 
0 
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where  0  and  1  are  subscripts  for  'up'  and  'down'  respec- 
tively and  W.(t)  is  the  renewal  density  for  the  event  "end 
of  stay  in  state  J".   Weiss  has  also  given  the  availability 
function  as: 


n(t)  =  R(t)  +  /   W,(t)  R(t-T)  dx 
0    1 


where  R(t)  is  the  reliability  function,  i.e., 


R(t)  =  /   fTT(x)  dx 
t 


B.   SOLUTION 

No  straight  forward  solution  of  these  equations  for 
n(t)  in  explicit  closed  form  has  been  found  except  in  the 
special  case  where  f„  and  rR  are  both  exponential  density 
functions.  However,  by  finding  the  Laplace  transforms  of 
the  probability  density  functions  and  the  renewal  density 
functions  the  Laplace  transform  of  the  availability  function 
was  found. 

1 .   Laplace  Transforms 

The  Laplace  transform  of  a  continuous  or  sectionally 
continuous  function  g(t),  t  in  [0,°°),  is 

«  OO  , 

g  (s)  =  /   e    g(t)dt  ;   s  is  a  complex  constant 
0 

whose  real  part  is  large  enough  that  the  integral  converges. 

The  Laplace  transform  of  f\j(t)  is 
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tt  00  j  00,  ■,    . 

fTT(s)  =  /   e"sr  fTT(t)dt  =  /   e"SL  A  e"AL  dt 
u      0        u        0 


=  tt-  /°°  U+s)  e"(A+s)t  dt  = 


A 


'*"•»/-         -v    x+s  ' 


:0t 


Similarly,  rn(s)  may  be  shown  to  equal 

The  Laplace  transforms  of  the  renewal  densities  are 

w*(s)  =  f*(s)  +  w*(s)  f*(s) 

w1(s)  =  wQ(s)  rD(s)  . 

These  imply 


»/  \  A(y+s)a 

wn(s)  = 


(A+s)(y+s)   -  Ay 


Noting  that  in  this  case 


wn(t) 
n(t)  = 


A 

the  Laplace  transform  of  the  availability  function  was 
found  to  be 


n*(s)  = (y+s)a 


(A+s)(y+s)a  -  Aya 


2 .   The  Inverse  Transform 

The  desired  solution  for  the  availability  function 

at 

is  the  inverse  transform  of  n  (s).   As  before,  no  method 
was  found  to  obtain  an  explicit  general  expression  for  n(t). 
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Therefore  a  numerical  inversion  procedure  was  used  to 
compute  values  of  n(t)  for  particular  values  of  the 
parameters  X,  u,  and  a  and  for  values  of  t=.05,  .10, 
.15j  •••*  1-0  mission  times. 

Several  values  were  chosen  for  the  parameters  X3 
u  and  a  over  a  range  that  was  judged  to  be  of  practical 
interest.   Those  values  and  the  corresponding  values  of 
ri(t)  are  shown  in  Appendix  A. 
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IV.   DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

The  computations  were  done  by  a  computer  program  written 
in  FORTRAN  IV  programming  language  and  run  on  the  IBM 
system/360,  model  67  computer  at  the  Naval  Postgraduate 
School  (see  program  listing,  Appendix  C).   The  program 
listing  contains  enough  explanatory  comments  to  enable  the 
reader  familiar  with  FORTRAN  to  follow  the  procedure. 

The  following  are  the  major  steps: 

1.  The  polynomial  in  the  denominator  of  ri  (s)  was 
expanded  to  the  form  aQ  +  a.s  +  a?s   +  ...  +  s     .   That  is, 

1  \  j  w    \0i    ,  a    „   /  ou   a-j  /,.    jy  \   J     ot+1 
(A+s)(y+s)   -  Ay   =  Z      (.)  y  d    (A+  a+\_< )    sJ  +  s 

J  —  ■* 
Note  that  an  equals  zero. 

2.  A  computer  library  subroutine  which  employs  the 
Newton-Raphson  successive  approximations  technique  was  used 
to  find  the  roots  of  this  polynomial. 

3.  After  the  roots  were  checked  to  see  that  they  were 
all  different  the  following  formula  was  used  in  computing 
the  values  of  n(t).   Often  called  Heaviside's  expansion 
formula  [4]  it  is  an  application  of  the  residue  theorem 
from  theory  of  complex  variables. 

Heaviside's  Expansion  Formula:   let  P(s)  and  Q(s)  be 
polynomials  where  P(s)  has  degree  less  than  that  of  Q(s). 

If  Q(s)  has  distinct  roots  A,  ,  k=l,2,...,n,  then  the 

P(s ) 
inverse  Laplace  transform  of  ^-7 — r  is 
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,-i/pooi     ?   p(Ak)     V 


By  application  of  this  formula 


L  (n  (s))  =  n(t)  =  s  ,+Ak+a(x+Ak)  e 
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APPENDIX  A 
RESULTS 

This  appendix  contains  tabulated  values  of  the  avail- 
ability function  for  the  indicated  sets  of  parameters  and 
the  indicated  time  in  mission-time  units. 

The  values  of  PZERO  given  are  such  that  PZERO  equals 
the  probability  of  a  repair  in  less  than  one  mission 
time ,  i.e., 

1    a      .,     . 

/  f    y-,  N  ;  ta  x   e  yt  dt  =  PZERO  . 
0   (a-1)! 

Zl  is  the  average  availability  over  the  first  mission 
time . 

zi  =  i  (n(0)^(1)  +  ?  n(.05i» 

dV)  d  i  =  l 

ZBAR  is  the  asymptotic  value  of  n(t)  for  large  t.   ZBAR 
is  also  the  long-run  average  availability. 


ZBAR  =     MTTF  1/X 


MTTP  +  MTR    1/X  +   a/u    y  +  aX 

Example:   For  a  system  whose  uptimes  have  the  expo- 
nential (A  =  .005)  distribution  and  whose  downtimes  have 
the  Erlang  (u  =  2.035  a  =2)  distribution,  the  true  system 
availability  at  a  half  mission  time  is  n ( . 5 )  "  .99777. 


17 


The  true  average  availability  over  the  first  mission 
time  is  Zl  =  .99792  and  the  long-run  average  availability 
is  ZBAR  =  .99510. 
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ALPHA=  2.0 

P7ERD=  0.60 

MU  = 

2.03 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.99751 

.99502 

.99007 

.93514 

.98023 

0.10 

.990^0 

.99901 

.99504 

.99011 

.93032 

.97063 

.96103 

0.15 

,99926 

.99852 

.99263 

.98531 

.97084 

.95653 

.94254 

Q.20 

.99902 

.99805 

.99027 

.Q8064 

. 96167 

.94307 

, 02434 

0.25 

.99379 

. 99759 

.98799 

.97614 

.95286 

.93016 

.90802 

0.30 

.99857 

.99714 

,98530 

.07130 

.94444 

.°1788 

. 89210 

0.35 

.99836 

.996  71 

.98369 

.96766 

.93642 

.90626 

. 37712 

0.40 

.99815 

.99631 

.98167 

.96370 

.92382 

.89530 

. 86308 

0.45 

.99795 

•  99591 

.97975 

,95995 

. o?164 

. 88501 

. 34997 

0.50 

.99777 

. 99554 

.97793 

.95639 

.91489 

.87533 

.33778 

0.55 

.99759 

.99519 

.97620 

.95303 

.90854 

.86640 

.32647 

0.60 

.99742 

.99485 

.97456 

.94987 

.90260 

.85304 

.81602 

0.65 

.99726 

.99454 

.97302 

.94689 

.89706 

.85029 

.80639 

0.70 

.99711 

.99424 

.°7157 

.94410 

.89189 

.84312 

.79754 

0.75 

.99697 

.99396 

.97021 

•  94149 

.88709 

.83650 

.78943 

o.ao 

.99684 

.99369 

.96894 

.93905 

.88264 

.330*0 

.73201 

0.85 

.99671 

.99344 

.96774 

.93678 

.87351 

. 8247° 

.77524 

0.90 

.99660 

.99321 

.96663 

,°3466 

.87470 

.81965 

,76908 

0.95 

.99649 

.99299 

.96559 

.93270 

.87119 

.31494 

.76349 

1.00 

.09639 

.99279 

.96462 

.93087 

.86795 

•  81064 

.75842 

Zl 


.99792   .99584   .97945   .95954   .92151   .88574   .85208 


ZB&R   .90510   .99Q24   .95305   .91031   .33539   .77136   .71731 
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ALPHA=  2.0 

PZER0=  0.70 

MU= 

2.45 

LAMBOA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.99  751 

.99502 

.99007 

.95515 

.98025 

0.10 

.99950 

.99901 

.99506 

.99014 

.98037 

.97071 

.96114 

0.15 

.99926 

.99853 

.99267 

.98539 

.97100 

.95682 

.94286 

0.20 

.99903 

.99806 

.99036 

.98082 

.96202 

.94360 

.92553 

0.25 

.99881 

.99762 

.98815 

.97646 

.95349 

.93109 

.90925 

0.30 

.99860 

.99719 

.98605 

.97231 

.94544 

.91935 

.89404 

0.35 

.°9839 

.99679 

.98406 

.96840 

.93788 

,90840 

.87992 

C.40 

.99820 

.99641 

.98219 

.96472 

.93081 

.89822 

,86689 

0.45 

.99802 

.99605 

.98043 

.96128 

.92424 

.88881 

.85492 

0.50 

.99785 

.99571 

.97878 

.95807 

.91816 

.88016 

.84397 

0.5  5 

.99770 

.99540 

.97725 

.95510 

.91254 

.87222 

.83400 

0.60 

.99755 

.99511 

.97582 

.95234 

.90738 

.86497 

.  824«>5 

0.65 

.99741 

.99484 

.97450 

.94979 

.90264 

.85836 

•  81676 

0.70 

.99729 

.99458 

.97328 

.94  745 

.89831 

.85236 

.80938 

0.75 

.99717 

.99435 

.97216 

.94529 

.89436 

.84693 

.80275 

0.80 

.99706 

•  99414 

.97112 

.94331 

.89076 

.84203 

.79681 

0.85 

.99696 

.99394 

.97017 

.94151 

.88750 

.83761 

.79150 

0.90 

.99687 

.99376 

.96930 

.93985 

,88454 

.83364 

. 78678 

0.95 

.99679 

.99359 

.96850 

.93835 

.88186 

.83008 

.78257 

1.00 

.99671 

.99343 

.96  777 

.93695 

.87945 

.82689 

.77885 

Zl 


.99803.   .99606   .98056   .96170   .92566   .89170   .85968 


ZBAR   .99593   .99190   .96078   .92453   .85965   .80328   .75385 
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ALPHA=  2.0 

P7ER0=  0.80 

MU  = 

3,00 

LAMBOA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

,400 

0.05 

.99975 

.99950 

.99751 

. 99503 

.99008 

.98516 

.98027 

0.10 

.99951 

.99901 

.99508 

.99018 

.98045 

.97083 

.96130 

0.15 

.9Q927 

.99854 

.99273 

.98552 

.97125 

.95719 

.94334 

0.20 

.99905 

. 99809 

.99050 

.98109 

.96255 

.94437 

.92656 

0.25 

.99883 

.99767 

.98839 

.97692 

.95441 

.93245 

.91103 

0.30 

.99863 

.99727 

.9864? 

.97304 

,94687 

.92146 

.89680 

0.35 

.99845 

.99690 

.98459 

.96944 

.93992 

.91141 

.88386 

0.40 

.99827 

.99655 

.98290 

.96613 

.93357 

.90226 

.87716 

0.45 

.99812 

.99624 

.98135 

.96310 

.92779 

.89400 

.86166 

0.50 

.99797 

.99595 

.97993 

.96034 

•9??56 

.88657 

.85227 

0.55 

.90784 

.99568 

.97863 

.95783 

.91784 

. 87991 

.84393 

0.60 

.99772 

.99544 

.97746 

.95556 

.91360 

.87398 

.83655 

0.65 

.99761 

.9952? 

.97640 

.95352 

.90981 

.86871 

.83004 

0.70 

.99751 

.99502 

.97544 

.95168 

.90643 

.86404 

.82433 

0.75 

.99742 

.99484 

.97458 

.95003 

.90342 

.85993 

.81933 

0.80 

.90734 

.99468 

.97381 

.94856 

.90076 

.85632 

.81498 

0.«5 

.99726 

.99454 

.97311 

.94724 

.89840 

.8531 5 

.81120 

0.90 

.99720 

.99441 

.97249 

.94607 

.89632 

.85037 

.8079? 

0.Q5 

.99714 

.99429 

.97194 

.94503 

.89449 

.84796 

.80509 

I. 00 

.99709 

.99419 

.97145 

.94411 

.89288 

.84586 

.80266 

Zl 


.99817   .99635   .98195   .96442   .93085   .89915   .86920 


ZBAR   .99668   .99338   .96774   .93750   .88235   .83333   .78947 


21 


ALPHA=  2.0 

PZFRO=  0.90 

MU= 

3.90 

LAMBDA 

TIMF 

.005 

•  010 

•  050 

.100 

.200 

.300 

•  400 

0.05 

,99975 

.99950 

•99752 

•99504 

.99011 

.98520 

.98031 

0.10 

.99951 

.99902 

.99512 

.99026 

.98061 

•97106 

.96161 

0.15 

.99928 

.99857 

.99285 

.98575 

.97171 

.95788 

.94426 

0.20 

.99907 

. 99814 

•99074 

. 98158 

.96352 

.94582 

.92846 

0.25 

.99888 

.99775 

.98881 

.97777 

.95608 

.93492 

,91428 

o.?o 

•99870 

.99740 

.98707 

.97433 

. 94°40 

.92520 

.90170 

0.35 

.9S854 

.99708 

.98550 

.97125 

.94346 

.91661 

.89066 

0.40 

.99840 

.99680 

.98410 

.96851 

.93823 

.90909 

.88106 

0.45 

.c9827 

. 99654 

.98287 

.96610 

.93364 

.90255 

. 87276 

0.50 

.99816 

.99632 

.98178 

.96399 

.92965 

.89689 

.86565 

0.55 

.99806 

.99613 

.98083 

.96215 

.92619 

.89204 

.85958 

0.60 

.99797 

•995°5 

.98000 

.06054 

.92321 

.88780 

.85444 

0.65 

.99790 

.99581 

.97928 

.95916 

.92066 

.88435 

.85011 

0.70 

.99783 

.99568 

.97865 

.95796 

.91847 

•  88136 

,84647 

0.75 

.99778 

.99556 

.97811 

.95694 

.91661 

.87884 

.84343 

0.80 

.99773 

.99547 

.97765 

.05605 

.91503 

•  87672 

.84090 

0.85 

.99769 

.99538 

.97725 

.95530 

.91370 

,87494 

.83831 

0.90 

.99765 

.99531 

.97691 

.95466 

.91257 

.87346 

.83708 

0.95 

.99762 

.99525 

•97662 

.95412 

.91162 

.87223 

.83566 

1.00 

.9S759 

.99520 

•97637 

.95365 

.91083 

.87121 

.83450 

Zl 

.99838 

.99676 

•98399 

.96341 

.93850 

.91013 

.88323 

78AR   .99744   .99490   .97500   .95122   .90698   .86667   .82979 


22 


ALPHA=  2.0 

PZERO=  0.95 

MU= 

4.80 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.9975? 

.99505 

.99013 

.98524 

.98037 

0.10 

•90O52 

.99903 

.99516 

.99035 

.98030 

.97134 

.96198 

0.15 

•Q9930 

.°9859 

.99998 

.98602 

.97225 

.95868 

.94531 

0.20 

.99010 

.99  8  20 

.99102 

.98212 

.96460 

.94741 

.93057 

C.25 

.90892 

.99784 

.98927 

.97868 

.95787 

.93757 

.91776 

0.30 

.90877 

.99753 

.98775 

.97567 

.95204 

.92910 

.90681 

0.35 

.9^863 

.99727 

.98643 

.97308 

.94705 

.92189 

.89755 

0.40 

.99P52 

.99704 

.98529 

,97086 

.94282 

.91581 

.88981 

0.45 

.99842 

.99684 

.98433 

.96899 

.93925 

.91074 

.88340 

0.50 

.09833 

.09667 

.98352 

.96741 

.93628 

.90654 

.87813 

0.55 

.99826 

.99653 

.98283 

.96603 

.93381 

.90308 

.87333 

0.60 

.99820 

.99641 

.98226 

.96498 

.93177 

.90025 

.87035 

0.65 

.99815 

.99631 

.98178 

.96407 

.93009 

.89795 

.86754 

0.70 

.00811 

.99623 

.98138 

.96331 

.92872 

.89609 

.86529 

0.75 

.99808 

.99616 

.98106 

.96269 

.92760 

.39458 

.86349 

0.80 

.99805 

.99611 

.98079 

.96218 

.92669 

.P9338 

.86207 

0.85 

.99303 

.99606 

.98056 

.96176 

.92596 

.89242 

.86095 

0.90 

.99801 

.99602 

.98038 

.96142 

.92537 

.89165 

.86007 

0.95 

.99799 

.99599 

.98023 

.96114 

.92490 

.89104 

.85938 

1.00 

.99798 

.99596 

.98011 

.96092 

.92452 

.89056 

.85885 

Zl    .99856   .99712   .93573   .97182   .94501   .91950   .89521 
ZBAR   .99792   .99585   .97959   .96000   .92308   .83889   .85714 


23 


Al_PHA  =  3.0 

P7ER0=  0.60 

MU= 

3.11 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

C.05 

.99975 

.09050 

.99750 

.99501 

.99005 

.98511 

.98020 

CIO 

.99950 

.99900 

,99502 

.99006 

.98022 

.97048 

.96083 

0.15 

.9Q925 

.99851 

.99255 

.98516 

.97054 

.95614 

.94195 

0.20 

.99Q01 

.99  802 

.99012 

.98034 

.96106 

.94217 

.92365 

0.25 

.Q9877 

.99753 

.98773 

.97562 

.95184 

. 92865 

.90603 

0.*0 

.99853 

.99706 

.98541 

.97103 

.94291 

.91563 

.88916 

0.35 

.Q9830 

.99661 

.98315 

.96659 

.93434 

.90320 

.87313 

0.40 

.99808 

.99616 

.98098 

.96233 

.92615 

.80138 

. 85799 

0.45 

.99787 

.99574 

.97889 

.95826 

.91837 

.88023 

. 84370 

0.50 

.99766 

.99534 

.97691 

.95439 

.91102 

.86978 

.83055 

0.55 

.09747 

•  99495 

.97503 

.95074 

.90413 

. 86002 

.81828 

0.60 

.99729 

.99459 

.97325 

.94730 

.89768 

.85097 

.80698 

0.65 

.90712 

.99424 

.97158 

.94408 

.89170 

. 84262 

.79664 

0.70 

,99695 

.99392 

.97002 

.94108 

.88616 

.83495 

.78721 

0.75 

.99680 

.99362 

.96857 

.93830 

.88106 

. 82795 

. 77868 

0.80 

.99666 

.99334 

.96722 

.93572 

.87638 

.82159 

.77099 

0.85 

.99653 

.993  0  8 

.96597 

.93334 

.87210 

.81583 

.76410 

0.90 

.99641 

.99284 

.96481 

.93116 

.86822 

.81064 

.75796 

0.95 

.99630 

.99261 

.96375 

.92917 

.86469 

.80599 

.75251 

1.00 

.99619 

.99241 

.96278 

.92735 

.86151 

.80183 

.74770 

Z\ 


.90782   .99564   .97849   .95767   .91797   .88071   .84572 


Z8AR   .99520   .99045   .95399   .91202   .83827   .77556   .72158 


2k 


ALPHA=  3.0 


P7FRO=  0.70 


Mt.l=  1,62 


LAMBDA 

TIMF 

.005 

.010 

.050 

.100 

.200 

.  300 

.  400 

0.05 

.9Q975 

.OQ950 

.09750 

.99501 

.990C5 

.98512 

. 93020 

0.10 

.99050 

. O9Q00 

.09502 

.09007 

.93023 

.97H49 

.96085 

0.15 

.99925 

.998  51 

.99256 

.98518 

. Q7059 

.95621 

.Q4205 

0.20 

.9Q901 

.99802 

.99015 

.98040 

.96120 

.04237 

.92392 

C.25 

.9^877 

.99755 

.98780 

.97576 

.95212 

.92907 

.90658 

0*30 

.9^854 

#qq?09 

.93554 

.S7120 

. 94343 

. 91640 

.80016 

0.35 

•cOR32 

.99665 

.983*6 

.06702 

.93517 

.00443 

.87474 

0.40 

.99811 

.99623 

.98130 

.96296 

.92739 

.8Q321 

. 860*8 

0.45 

.CQ791 

.Q0583 

.97934 

.95915 

.92011 

•  88279 

,3471? 

0.50 

.99773 

.99546 

.07751 

.95559 

.°1335 

.87318 

. 834Q8 

0.55 

,09755 

.99511 

.97581 

.95227 

.90711 

. 364*n 

•  82304 

0.60 

.09739 

.99478 

.07422 

.94921 

. 00 140 

.85637 

.81397 

0.65 

.09724 

•  99448 

.97277 

.94640 

.89619 

.84913 

.8050* 

0.70 

.09710 

.99420 

.97143 

.94384 

.80147 

. 84263 

.70708 

0.75 

•9Q697 

.99395 

.97020 

.94150 

.83721 

.83682 

.790Q4 

0.80 

.9Q685 

.99372 

.96909 

.93939 

.8^340 

.83167 

. 78386 

0.85 

.99675 

.09351 

.06808 

.93748 

, 87909 

.82711 

,77»46 

0.90 

.99665 

.99332 

.96717 

.93577 

.87696 

.82311 

.77*77 

0.95 

.90656 

.99314 

.96635 

.93423 

. 37429 

.81962 

.76973 

1.00 

.9C649 

•  99299 

.06562 

.93287 

.3719* 

.81658 

.76627 

71    .Q0791   .99583   .97940   .95945   .92138   .38562   .35200 
7BAR   .99587   .99178   .96021   .92347   .85732   .80088   .75104 


25 


ALPHA*  3.0 

<>ZER0=  0.80 

MU= 

4.30 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.99750 

.99501 

.99005 

.98512 

.98021 

0.10 

.99950 

.99900 

.99503 

.99008 

.98025 

.97052 

.96089 

0.15 

.99926 

.<?9951 

.99259 

.98523 

.97067 

.95634 

.94221 

0,20 

.99902 

.99803 

.99021 

.98052 

.96142 

.94270 

.92436 

0.25 

.99879 

.99757 

.98792 

.97600 

.95259 

.92976 

.90750 

0.3  0 

.99857 

.99713 

.98575 

.97171 

.94426 

.91762 

.89177 

0.35 

.99836 

.99672 

.  ?8370 

.96768 

.93648 

.90636 

.87729 

0.40 

•  99tU6 

.99633 

.98179 

.96394 

.92930 

.89603 

.86407 

C.45 

.99798 

.99597 

.98002 

.96049 

.92274 

.88665 

.85215 

0.50 

.99782 

.99564 

.97841 

.95735 

.91678 

.87820 

•84151 

0.55 

.95766 

.99534 

.97694 

.95450 

.91143 

.87068 

.83210 

0.60 

.99753 

•  ^  9506 

.97561 

.95193 

.90666 

.86402 

.82385 

0.65 

.99740 

.99482 

.97441 

.94964 

.90243 

.85818 

. 81668 

0.70 

.99729 

.99460 

.97335 

.94760 

.89872 

.85309 

.81050 

0.75 

.99719 

.99440 

.97240 

.94581 

.89547 

.84869 

.80522 

0.80 

.99711 

.99422 

.97157 

.94423 

.89264 

.84491 

.80074 

0.85 

.99703 

.99407 

.97083 

.94284 

.89020 

.84169 

.79697 

0.90 

.99696 

.99393 

.97019 

.94164 

.88810 

.83996 

.79383 

0.95 

.99690 

.99381 

.96963 

.94060 

.88631 

.83666 

.79122 

1.00 

.99685 

.99371 

.96914 

.93970 

.88479 

.83474 

.78908 

Zl    .99803   .99608   .98062   .96183   .92595   .89218   .86038 
ZBAR   .99652   .99307   .96629   .93478   .87755   .82692   .78182 


26 


ALPHA=  3.0 

PZERO=  0.90 

M(J= 

5.35 

LAMBDA 

TIMF 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

• 00050 

.99750 

.99502 

.99006 

.98512 

.98021 

0.10 

.9Q950 

.99901 

.99504 

,99010 

.98029 

.97Q58 

.96097 

0.15 

.99926 

.9Q852 

.99263 

.98  531 

. 97084 

, 95659 

.94255 

0.20 

.9Q903 

.99  806 

„ qQ03? 

.98074 

. 96187 

.94337 

.92524 

0.25 

.99881 

.99762 

.98815 

.07646 

.95349 

.93110 

.  90927 

0.30 

.cq861 

.99721 

.98614 

. 97250 

.94581 

.91992 

.89479 

0.  35 

.99842 

. 99684 

.98431 

.96890 

.93887 

.90988 

.88189 

0.40 

.99825 

.Q9650 

.98  26  6 

.96566 

,93268 

. 90100 

.87055 

C.45 

.99810 

. 99620 

.98119 

.96280 

.92^24 

.89324 

.86074 

0.50 

.C9797 

.99594 

.97990 

.96029 

.92251 

. 88656 

. 85236 

0.  55 

.99785 

.99571 

,97877 

. 95810 

.91843 

.83087 

•  84528 

0.50 

.99775 

.09551 

.97779 

.95623 

.91496 

.87606 

.83938 

0.65 

.9Q766 

.99533 

.97695 

.95462 

.91203 

.87205 

. P3450 

0.70 

.99759 

.99518 

.97623 

.95326 

.909*7 

.86873 

.83052 

0.75 

.Q9752 

.99506 

.97563 

.95211 

.90753 

.86601 

.82731 

C.PO 

.99747 

.99495 

.97512 

.95115 

.90585 

. 86379 

.82474 

0.85 

.Q9742 

.99486 

.97469 

.95036 

.90447 

.86201 

.82^71 

0.90 

.99^39 

.99478 

.97433 

.94970 

.90334 

.86059 

.32112 

0.Q5 

.99735 

. 99472 

.97403 

. 94915 

.90244 

.85947 

.81900 

1.00 

.99733 

.99466 

.97379 

.94871 

.90171 

.85858 

,81896 

Zl    .99822   .99644   .98241   .965^4   .93266   .901*1   .87268 
ZBAR   .99720   .09442   .97273   ,94690   .8Q916   . 8560C   .81679 


27 


ALPHA=  3.0 

PZER0=  0.95 

MU= 

6.30 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.99751 

.99502 

.99006 

.98513 

.98022 

0.10 

.99950 

.99901 

.99505 

.99012 

.98034 

.97066 

.96107 

0.15 

.99927 

.99853 

.99268 

.98541 

.97105 

.95689 

.94295 

0.20 

.09904 

.99  80  8 

.99045 

.98100 

.96238 

. 94413 

.92624 

0.25 

.99883 

.99767 

.98841 

.97696 

.95448 

•93?57 

.91119 

0.30 

.99865 

.99730 

.98656 

.97333 

.94  744 

.92233 

. 89796 

0.35 

.99848 

•  99697 

.98493 

.97013 

.94129 

.91345 

.88656 

0.40 

.99834 

.99668 

.98352 

.96736 

.93600 

.90537 

.87692 

0.^5 

.99821 

.99643 

.98230 

.96500 

.93153 

. 89952 

.86891 

C.50 

.99811 

.99622 

.98127 

.96300 

.92779 

.89426 

,86235 

0.55 

.99802 

.99604 

.98041 

.96134 

.92470 

.88997 

.85705 

0.60 

.99794 

.99589 

.97969 

.95997 

.92219 

.88652 

.85284 

0.65 

.99788 

.99577 

.97910 

.95884 

.92015 

.88377 

.84954 

0.70 

.99783 

.99567 

.97862 

.95793 

.91853 

.88160 

.84698 

0.75 

.99779 

.99559 

.97823 

.95720 

.91724 

.87991 

.84502 

0.80 

.99776 

.99552 

.97792 

.95661 

.91623 

.87861 

.84354 

0.85 

.99773 

.99547 

.97766 

.95615 

.91544 

.87762 

.84245 

0.90 

.99771 

.99542 

.97746 

.95578 

.91484 

.87688 

.84164 

0.95 

.99769 

.^9539 

.97730 

.95549 

.91437 

.87632 

.84106 

1.00 

.99768 

.99536 

.97718 

.95527 

.91402 

.87591 

.84065 

11    .99837   .99674   .98388   .96821   .93815   .90970   .88274 
7BAR   .99762   .99526   .97674   .95455   .91304   .87500   .84000 


28 


ALPHA=  4.0 

PZERO=  0.60 

MU  = 

4.20 

LAMBDA 

TIMF 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.QQ975 

.99950 

.99750 

.99501 

.99005 

.98511 

.98020 

0.10 

.99950 

# QQ90  0 

,99501 

.09005 

.93020 

.97045 

.96080 

0.1  5 

.OQQ?5 

. 9985^ 

.99253 

.98512 

.97047 

.95604 

.94182 

0.  2  0 

.99900 

.99801 

.99007 

.98025 

.9608P 

.94191 

.92330 

0.?5 

.99876 

.907  5  2 

.98764 

.97544 

, 9514« 

.92812 

.90533 

0.  3<^ 

. 9CR5? 

.99703 

.98525 

.97072 

.04231 

.9147^ 

. 887DQ 

0.35 

.90823 

.09656 

.98292 

.96613 

.93343 

.90135 

.87137 

o.^p 

.09805 

•  99610 

.98065 

.96168 

.92488 

.88052 

,«5556 

0.45 

.09732 

.99565 

.97847 

.Q5741 

.9167! 

.87781 

.84064 

0,50 

.99761 

.09523 

.97637 

.95334 

.90397 

. 8667* 

. 32667 

0.55 

.99741 

.99482 

.97430 

.94948 

.00163 

.85647 

.81371 

0.60 

•°9721 

.99444 

.97251 

.94584 

.89488 

.84601 

.*C177 

0.65 

.OC703 

.99^07 

.97075 

.94245 

.88856 

.83811 

.79087 

0.70 

.99686 

.99373 

.96911 

,939?9 

.8*274 

.83007 

.78101 

0.75 

.Q°670 

.9934? 

.96753 

.93638 

.87742 

.«2278 

.77214 

o.ao 

.99656 

.99313 

.9661 8 

.93370 

.87258 

.81622 

.76425 

0.85 

.90642 

.09286 

.96480 

,93126 

.86  821 

.81036 

.75727 

0.90 

.00630 

.99261 

.96371 

.02905 

. 86428 

.80515 

.75117 

0.9  5 

.99618 

.99238 

.96264 

.92  704 

.86078 

.80057 

.74536 

1.00 

.9Q608 

.99218 

.96168 

.92524 

.85767 

.70657 

.74130 

71 


.99776   .99553   .97795   .95661   .91597   .87736   .84212 


7BAR   ,99526   .99057   .95455   .91304   .84000   .77778   .72414 
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ALPHA=  4.0 

PZFR0=  0.70 

MU  = 

4.76 

LAMBDA 

TIME 

.005 

.010 

.050 

.  100 

.200 

.300 

.400 

0.05 

.99975 

•99O50 

.99750 

.99501 

.99005 

.98511 

.98020 

0.10 

.99950 

.99900 

.99501 

.99005 

.98020 

.97045 

.96080 

0.15 

.Q9925 

.99850 

.99254 

.98513 

.97049 

.95606 

.94184 

0.20 

.99900 

.99801 

.99009 

.98027 

.96093 

.94198 

.92340 

0.25 

.99876 

.90752 

.98767 

.97550 

.95161 

.92831 

.90553 

0.30 

.99852 

.99705 

.98532 

.97086 

.94258 

.91513 

.88851 

0.35 

.99829 

.99658 

.98304 

.96637 

.93390 

.90256 

.87229 

0.40 

.99807 

.99614 

.98085 

.96208 

.92565 

.89066 

.85705 

0.45 

.99785 

.99571 

.97876 

.95800 

.91787 

.87952 

.84287 

0.50 

.99765 

.99531 

.97680 

.^5418 

.91061 

.86919 

.32982 

0.55 

.99746 

.99494 

.97496 

.95061 

.90389 

.85971 

.81791 

0.60 

.99729 

. 99459 

.97325 

.94731 

.89773 

.85108 

.80718 

0.65 

.99713 

.99426 

.97168 

.94429 

.89213 

.84330 

.79759 

0.7Q 

.99698 

.99396 

.97025 

.94153 

.88708 

.83636 

.78911 

0.75 

.99684 

.99369 

.96894 

.93904 

.88256 

.83021 

.78169 

0.80 

.99672 

.99345 

.96776 

.93681 

.37855 

.82481 

.77526 

0.85 

.99661 

.99323 

.96671 

.93482 

.87501 

.82012 

.76974 

0.90 

.99651 

.99303 

.96576 

.93305 

.87192 

.81608 

.76506 

0.95 

.99642 

.99285 

.96493 

.93150 

. 86924 

.81263 

.76113 

1.00 

.99634 

. 99269 

.96419 

.93013 

.86692 

.80970 

0 

.75787 

Zl 

.99784 

.99568 

.97870 

.95807 

.91377 

.88191 

. 84730 

ZBAR   .99582   .99167   .95968   .92248   .85612   .79866   .74843 
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ALPHA=  4.0 


p/FRH=  0.B0 


Ml)=  5.57 


LAMBDA 

TIME 

.005 

.010 

.05^ 

.100 

.200 

.300 

.400 

0.0  5 

.99975 

#  09050 

.99750 

.99501 

.99005 

.98511 

. 98020 

0.  10 

.09950 

. O9900 

.09502 

.99006 

.98021 

.97046 

,06081 

0,1  5 

.90075 

•  09  P  5  0 

.99*55 

.08515 

.97051 

. 95610 

, 9410Q 

O.*0 

.99901 

•c9801 

.90011 

.08032 

.96103 

. 94212 

.92359 

0.25 

.°9877 

.09753 

.98773 

.07562 

.95184 

.07865 

.00603 

0.30 

•°9853 

.99707 

.08544 

. 97109 

. 0^304 

.°1582 

,Q«041 

0.3  5 

.99831 

•  09663 

.98324 

.9667  8 

.93471 

.90^76 

.37337 

0.40 

.99810 

.09620 

.98118 

.96273 

.92694 

.30756 

.85954 

0.45 

.99790 

•  90581 

.07025 

.958Q7 

,91976 

, 88220 

.S4650 

0.50 

.99772 

. 99545 

.97747 

.95551 

.91321 

.87300 

.83478 

0.55 

.99755 

.90512 

,°7585 

.95236 

.90  731 

.8  647  0 

.82440 

0.60 

.90740 

.90481 

.97438 

.04953 

.Q0205 

.85736 

.815*2 

0.65 

.09727 

. 09454 

.07307 

.94701 

.89740 

.35096 

. 80748 

0.70 

.99714 

.99430 

.07190 

.94478 

.P0334 

. 84543 

. 80078 

0.75 

.Q9704 

. 99409 

.97087 

.94283 

.88983 

.84070 

.79515 

0.80 

.Q9694 

.99390 

.06007 

.94113 

. 8*682 

. 33671 

. 79046 

0.85 

.99686 

. 9937* 

,96910 

.93967 

.88427 

.83338 

."'9661 

0.90 

.99679 

.993  59 

,96852 

.93842 

.3Q212 

.33062 

.78*50 

0.95 

.09673 

.99347 

.96794 

.03735 

. 93032 

. 32837 

.78102 

1.00 

.99667 

.99^36 

.96745 

,93645 

.878«4 

.32656 

.77908 

Z! 


.99795   .99590   .97974   .96013   .92271   .88757   .85454 


ZBA9   .99639   .99281   .965^3   .03243   .87342   .82143   .77523 
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ALPHA=  4.0 

PZER0=  0.90 

MU= 

6.70 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

.99950 

.99750 

.99501 

.99005 

.98511 

.98020 

0.10 

.99950 

.99900 

.99502 

.99006 

.98022 

.97048 

.96083 

0.15 

.99925 

.99851 

.99256 

.98518 

.97058 

.95619 

.94202 

0.20 

.99901 

.99802 

.99016 

.98042 

.96124 

.94243 

.92400 

0.25 

.99878 

.99756 

.98786 

.97587 

.95233 

.92  93  8 

.90700 

0.30 

.99856 

.99712 

.98568 

.97157 

.94400 

.91724 

.89128 

0.35 

.99835 

.99671 

.98366 

.96760 

.93633 

.90615 

.87701 

0.40 

.99816 

. 99633 

.98181 

.96399 

.92940 

.89620 

.86430 

0.45 

.99799 

.99599 

.98015 

.96075 

.92325 

. 88743 

.85320 

0.50 

.99784 

.99569 

.97863 

.95789 

.91787 

.87983 

.84367 

0.5  5 

.99771 

.99543 

.97739 

.95541 

.91323 

.87335 

.83563 

0.60 

.99760 

.99520 

.97628 

.95327 

.90930 

. 86791 

.82896 

0.65 

.99750 

.99500 

.97533 

.95146 

.90599 

.86341 

.82353 

0.70 

.99741 

.99483 

.97453 

.94993 

.90326 

.85975 

.81917 

0.75 

.99734 

• $9469 

.97385 

.94867 

.90103 

.85680 

.81574 

0.80 

.99728 

.99457 

.97329 

.94763 

.89922 

.85448 

.81308 

0.85 

.99723 

.99448 

.97283 

.94678 

.89778 

.85266 

.81107 

0.90 

.9<5719 

. 99440 

.97246 

. 94609 

.89665 

.85127 

.80958 

0.95 

.99716 

.99433 

.97215 

.94554 

.89576 

.85022 

.80850 

1.00 

.99713 

.99428 

.97191 

.94510 

.89508 

. 84944 

.80775 

Zl    .99811   .99623   .98136   .96328   .92375   .89625   .86563 
7BAR   .99702   .99407   .97101   .94366   .89333   .84810   .80723 
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Al_PHA  =  4.0 

P7ER0=  0.95 

MU  = 

7.80 

LAMBDA 

TIME 

.005 

.010 

.050 

.100 

.200 

.300 

.400 

0.05 

.99975 

. 99950 

.99750 

.99501 

.99005 

.98511 

.98020 

0.10 

.99950 

.99900 

.a9502 

.99007 

. 98024 

.97050 

.96086 

0.15 

.9Q926 

.99851 

.99258 

.98522 

.Q7066 

.95632 

.94219 

0.20 

.99902 

• 99304 

.99023 

.98056 

.°6150 

.9428^ 

.92453 

0.25 

,o°879 

.997*9 

.98801 

.97618 

.95294 

. 93029 

.90820 

0.  30 

.°9859 

.99718 

.98597 

.97215 

.9451? 

.91892 

.89350 

0.35 

.99P40 

.99680 

.98413 

.96854 

.93817 

.90887 

. 88058 

0.40 

.99823 

.99647 

.98251 

.96537 

. 93211 

.90018 

.86952 

0.^5 

. QQ809 

. 99619 

.98110 

.96263 

.92694 

.89233 

.86025 

0.50 

.99797 

.09*94 

.97991 

.96032 

•  92260 

. 8«675 

.85267 

0.  5  5 

.99786 

. 99573 

.97891 

.95840 

.91904 

.83181 

,34659 

0.60 

.99778 

.99556 

.97808 

.95682 

.91616 

. 87738 

.84182 

0.65 

.99771 

.9954? 

.97741 

.95554 

.91337 

.87480 

.83816 

0.70 

.99765 

.99531 

.97687 

.95452 

.91207 

.87244 

.33541 

0.75 

.99760 

.99522 

.97644 

.95372 

.91  069 

. 87065 

. 83339 

0.30 

.90757 

.99515 

.97610 

.95309 

.90963 

.36934 

.83195 

0.85 

.99754 

.99509 

.97583 

.95261 

.90884 

.86838 

.83095 

0.90 

.99752 

.99504 

.97563 

.95224 

.90  8  26 

. 86771 

.83023 

0.95 

.99750 

.99501 

.97547 

.95196 

.90783 

.86724 

•«2935 

1.00 

.99749 

.99498 

.97535 

.95175 

.90753 

. 86693 

.82960 

Zl    .99825   .99651   .98277   .96604   .93403   .90382   .87528 
ZBAR   .99744   .99490   .97500   .95122   .90698   .86667   .32979 
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APPENDIX  B 

EXAMPLE  PROBLEM 

Suppose  a  system  consisting  of  15  components  is  connected 
in  series.   Each  component,  operating  independently  of  the 
others,  has  an  exponential  failure  density  with  parameter  A 
and  a  special  Erlang  repair  density  with  parameters  u  and 
a.   Assume  that  uptimes  and  downtimes  are  all  independent. 

a.  What  is  the  average  availability  of  the  system  over 
the  first  mission  time  if; 

(1)  A  =  .005  failures  per  mission  time,  u  =  4.8 
repairs  per  mission  time,  and  a  =  2? 

(2)  A  =  .40  failures  per  mission  time,  u  =  4.2 
repairs  per  mission  time,  and  a  =  4? 

b.  What  is  the  long-run  average  availability  of  the 
system  if; 

(1)  the  parameters  are  the  same  as  in  a(l)? 

(2)  the  parameters  are  the  same  as  in  a(2)? 
Solution: 

a. 

(1)   The  average  availability  of  a  single  com- 
ponent over  the  first  mission  time  is  Zl  =  .99856  (from 
Appendix  A).   Letting  T\      be  the  average  system  availability 
over  the  first  mission  time; 


ns 


(Zl)15  =  (.99856)15  =  .9786  . 
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(2)    In   this    case    Zl  =  .8^212    and   n     =  (.84212)15  =   .0760 

s  

b. 

(1)  The  long-run  average  availability  of  a  single 
component  is  ZBAR  =  .99792. 

Letting  n   be  the  long-run  average  availability  of  the 
system, 

if  =  (ZBAR)15  =  (.99792)15  =  .9693  • 
s 

This  compares  with  the  average  system  availability  over 
the  first  mission  time  of  .9786. 

(2)  In  this  case  ZBAR  =  .  724l4  and  r\      =  .0079  . 
This  compares  with  the  average  system  availability  over 

the  first  mission  time  of  .076. 

The  above  comparisons  indicate  the  use  of  the  long-run 
average  availability  to  approximate  the  system  availability 
over  the  first  mission  time  is  a  conservative  practice. 
These  cases  and  others  also  indicate  that  the  lower  the 
value  of  ZBAR  the  more  conservative  is  the  practice. 
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APPENDIX  C 


FORTRAN  IV  COMPUTER  PROGRAM  FOR  THE 
INVERSION  OF  n*(s)  BY  NUMERICAL  METHODS 


C 
C 

c 
c 
c 

c 

c 


c 
c 


THIS    PROGRAM    COMPUTES    VALUES    OF     SYSTEM    AVAILABILITY    BY 
INVERTING    THE    LAPLACE    TRANSFORM    OF    THE    AVAILABILITY 
FUNCTION, 

IMPLICIT    REAL*8(  A-  H),  REAL*  8(0-7.  ) 
COMPLEX* 16    S, ETA, EX  PON, CONST, TOP, BOTTOM 
DIMENSION    XCOFC201 ,COF( 20) , R0HTR(2  0) ,RHGTI ( 20 ) , S( 20 ) , 
1EXPON(20),RCTA(20,7) ,YLAM0A(7>, ZBAP(7),Z1<7) ,ETA(  ?0) 
DATA    YL AMD A /O. 0^5, 0.01 0,0, 05, 0.1 ,0.2,0,3,0.4/ 
JCASE=0 


1 

10 


ICASF    IS    THF    NUMBER 
PEAD(5,1  UCASE 
FORMAT( 12) 
READC5,10)ALPHA,XMU,PZER0 

F0RMAT(F3.0,2F10.3) 
JCASE=JCASF+1 

WRITE(6,6) JCASE 

6  FORMAT* 1H1 ,45X,»CASE« ,14) 
DO  800  K=l,7 
XLAMDA=YLAMDA(K) 

WRITE(6,7) 

7  F0RMAT(//10X, • INPUT    DATA') 
WRITE(6,8) ALPHA,     XMU,XLAMDA 

8  FORMATI/10X, 'ALPHA    =     »,F3,0,5X,« 
l'LAMBDA    =• ,F10.3) 

PUT    ZEROS     IN    ALL    ARRAYS. 
DO    80    L=l,20 
XCOF(L )=0,0 
COF(L)=0,0 
RO0TR(L)=0.0 
ROOTKL  )=0,0 
S  (L)  =  (0. 0,0,0) 
ETA(L)=(0. 0,0.0) 
EXPON(L )=(0. 0,0.0) 


HF    CASES    TO    RE    PROCFSSEO    THIS    RUN. 


MU    =• ,F10. 3,5X, 


80 


M=ALPHA 

N  =  M+1 

IF(M,NE.ALPHA)GO 


TO    40 


COMPUTE  COEFFICIENTS  OF  THE  POLYNOMIAL, 
CALL  COEFFI ALPHA, XMU, XL  AMD A, XCOF ) 

FIND  THE  ROOTS  OF  THE  POLYNOMIAL, 
OPOLRT  IS  A  PROGRAM  LIBRARY  SUBROUTINE. 

CALL    DP0LRT(XC0F,C0F,N,R0OTR,RO0TI , IER ) 


60    WRITE(6,6i ) IER 
IER    IS    AN    ERROR    CODE 


61 

18 

19 
14 

16 


FROM  SUBROUTINE  DEPOLRT, 


IER=»,I2) 


OF    THE    POLYNOMIAL. • I 


FORMAT(//10X,  • 

WRITE(6,18) 

F0RMAT(/13X, 'ROOTS 

WRITE(6,19  ) 

FORMAT( /15X,«REAL» 

WRITE(6,16)  (ROOTRI  I  ),ROOTI( I ), 

FORMAT(10X,2F12.5) 


6X 


IMAGINARY*/     ) 
1=1, N) 


CHFCK  FOR  MULTIPLE  RO^TS. 

CALL    CKMULT<ROOTR,ROOTI,M) 
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c 

C   PUT  THE  ROOTS  INTO  THE  COMPLEX  ARRAY  S(20). 
DO  100  1=1, N 
100  S( I ) =0CMPLX< ROOTR  ( I  ) , ROOT  I ( I) ) 
C 
f   BFGIN  INVERSION  BY  HFAVISIDE'S  EXPANSION  FORMULA. 

c 

DO  150  1=1, N 

TOP=XMU*S( I ) 

R0TT0M=Tnp4-V*  (XLAMDA+S<  I  )  ) 

CONST=TOP/RPTTOM 

DO  150  J=l,20 

EXPON( J)=CONST*CDEXP(  J*S( I )*0.05l 
150  ETA( J)=ETA< J)+EXPON(J) 
C 

WRITE(6,160) 
1*0    F0RMAT(///10X, 'VALUES    OF    THE    AVAILABILITY    FUNCTION    AT' 
1'     .05    MISSION-TIME     INCREMENTS. •/ » 

DO    170    1=1 , ?0 

T=0.05*I 

170  WRITE(6,171)T,FTA( I > 

171  FORMAT! 10X,F4. 2,2F10. 6) 
C 

SUM=0.0 
DO    155    L=1,?0 
RETA(L,K)  =  PEAL(FTA(L)  ) 
155    SUM=SUM+RETA(L,K) 

T0TSUM=SUM*SUM+1,-RETA(20,K) 

C 

C   71  IS  THE  AVERAGE  AVAILABILITY  HVER  THE  FIRST  MISSION-TIM 

Zl(K)=T0TSUM/40. 
C 
C   ZRAR  IS  THE  LONG  RUN  TIME  AVERAGE  AVAILABILITY. 

300    ZBAR(K  )=XMU/<  ALPHA*XLA MDA+XMU) 
C 

CALL  OUTPUT (YL AMD A, ALPHA, P ZERO, XMU,Z1, 7  BAR, RET  A) 
C 

30  IF(JCASE-ICASE)5,50,*0 

40  WRITE(6,41 ) 

41  FORMATdOX  , 'ERROR-  ALPHA  MUST  BE  AN  INTEGER'  I 
GO  TO  30 

50  STOP 
END 

SUBROUTINE  CKMULT ( ROOTR , ROOT  I , M ) 

IMPLICIT  REAL*8( A-H) ,REAL*8(0-Z) 

DIMENSION  S  (1),  ROOTR  (1  I.ROOTK  1  ) 
C 
C   THIS  SUBROUTINE  CHECKS  FOR  MULTIPLICITIES  AMONG  THE  ROOTS 


C 


MULT=1 

DO  50  1=1, M 

K=M+l-t 

DO  50  J=1,K 

50  I  F(  ROOTR  (I  ).FQ.ROOTR(I+J  I.  AND- ROOT  I  (  I  J.EQ.ROOTI  (  I  ♦  J  I  ) 
1MULT=MULT+1 

IF(MULT.NE.l)  GO  TO  60 

WRITE(6,51  ) 

51  FORMAT(///10X, 'THERE  ARE  NO  MULTIPLE  ROOTS.') 
RETURN 

6C  WRITE(6,61) 

61  F0RMAT(///10X, 'THE  POLYNOMIAL  HAS  MULTIPLE  ROOTS.  THIS 
l'PROGRAM  WILL  NOT  ACCURATELY  HANDLE  THIS  CASE.') 

RETURN 

END 
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SUBROUTINE  COEFF ( ALPHA , XMU, XL AMD  A, XCOE ) 

IMPLICIT    RE/>L*8<  A-H) , R EAL*8 ( 0- Z  ) 

DIMENSION    XCOF(l) 
C 

C       THIS    SUBROUTINF    COMPUTES    THE    COEFFICIENTS    OF     THE    EXPANDED 
C       POLYNOMIAL     IN    THE    DENOMINATOR    OF    THE    LAPLACE    TRANSFORM. 
C 

M=ALPHA 

N2=M+2 

X=ALPHA+1.0 

XC0F(1)=0.0 

XCOF(M4-2)=1.0 

CHOOZ=M 

DO    20    1=1, M 

BKT=XLAMDA+(XMU*I ) /(X-I ) 

XC0F(H-1)  =  CH00Z*BKT*XMU**(M-I) 
20    CHOOZ=CHOOZ*(M-I )/(I+l ) 
C 

C   THF  COEFFICIENTS  ARE  NOW  IN  THE  ARRAY, XCOF. 
C 

WRITE(6,17) 
17  FORMAT(//10X, 'COEFFICIENTS  OF  THE  POLYNOMIAL1  I 

WRITE (6, 12)  (XCOF (I )  ,I  =  1,N2) 
12  FORMAT!  /10X, 10F12.6) 

RETURN 

END 

SUBROUTINE  OUTPUT (YLAMDA, ALPHA , PZE RO, XMU, Z 1 , ZBAR ,R ETA ) 
IMPLICIT  REAL*3(A-H),REAL*8(0-Z) 
DIMENSION  YLAMDA(1),Z1( 1» , ZBAR ( 1 ) , RFTA ( 20, 7 ) 
WRITEC 6,1) ALPHA, PZERO, XMU 
1  FORMAT ( 1H1 ,////////25X,» ALPHA=« ,F4,1 ,5X,»PZER0=» , F5.2, 
15X,«MU=,,F5.2) 
WRITE(6,11) 

11  FORMATC //46X, •LAMBDA1) 
WRITE (6, 2) (YLAMDA(L),L=1,7) 

7    FORMAT! /22X ,F4. 3 ,6 ( 4X , F4.3 ) ) 
WRITE(6,12) 

12  F0RMAT(15X, 'TIME'/ ) 
DO  50  1=1,20 
T=0.05*I 

50  WRITE(6,3)T,(RETA(I, J) , J=l,7) 

3  FORMAT! 15X,F4.2,7(2X,F6. 5) /) 
C 

WRITE(6,A) (Zl( I), 1=1, 7) 

4  FORMAT (/16X,»Z1» ,1X,7(2X,F6.5)//) 
WRITE (6, 5)  ( ZBAR(  I) ,  1  =  1,7) 

5  F0RMAT(15X,«ZBAR« ,7(2X,F6.5) ) 
RETURN 

END 
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